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Abstract 

The coexistence of superconducting and charge-density-wave order in the half-filled 
attractive Hubbard model is interpreted as a consequence of the pseudospin SU(2) 
symmetry spontaneously broken to a 'hidden' subgroup U(l). By topological argu- 
ments we show that there must exist new skyrmion textures associated with this 
symmetry breakdown. This fact is illustrated via a non-linear a-model. Unlike the 
spin textures previously known in an antiferromagnetic background, doping the model 
away from half-filling leads the new skyrmions to unwrap. 
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Recently, the Hubbard model has been a subject of intense study for its possible 
applications to high T c superconductivity [1]. It represents the simplest theoretical 
framework for describing strongly correlated electron systems. Despite its simple 
appearance, it is however quite difficult to analyze beyond one dimension. The elim- 
ination of phonon or other internal degrees of freedom can give rise to an effective 
attraction between the remaining electrons [2]. In this paper, we shall restrict our- 
selves to the Hubbard model of attractive interaction, which has been frequently 
applied to such compounds as bismuth oxide superconductors [3] and many others 
[4]. Numerical studies [5] showed that for dimension d > 2 it has a ground state of 
both charge-density wave (CDW) order and superconductivity (SC) when half-filled. 
We shall show how an effective action for low energy excitations around such a ground 
state can be derived directly from the spontaneous breakdown of the SO (4) symme- 
try of the model [6,7,8], and that there must exist a new kind of skyrmion textures 
associated with this breakdown. 

Let us consider the ci-dimensional (d = 2 or 3) Hubbard model on a square or cubic 
lattice. The Lagrangian reads (summation convention on repeated indices assumed) 

L = ^2cl(x,t)id t c a (x,t) +t ^ [4(x,t)c T (x',t) + /i.c.] -U^n t (x,i)n|(x,t) , 

x {xx'J x 

(1) 

where c^x, t) and c£.(x, t) respectively annihilates and creates electrons with spin a 
(=|, |) at site x = (x±, ■ ■ ■ ,Xd) and time t, and n CT (x) = 4(x)c (T (x) (no summation). 
Hopping is restricted to nearest neighbors, as indicated by the bracket (xx'), with t 
the (constant) transferring matrix element. The constant U is the on-site Hubbard 
interaction. This is a standard theory for electrons in well- localized atomic orbitals 
with a probability t for transitions between neighboring atoms. By now it is well 
known that the Hubbard model at half-filling has SO(4)~ SU(2) c x SV(2) S /Z 2 sym- 
metry [6,7,8], where SU(2)s corresponds to spin rotational invariance and SU(2)c 
(sometimes called pseudospin group [9]) contains the ordinary U(l)c charge symme- 
try as a subgroup. To see the SO(4) symmetry of the Lagrangian (1) in a manifest 
way, we use the matrix representation as introduced by Affleck [6] and Schulz [8], 

1 ] \ (-) x c|(x,t) -R x c{(x,t) ) ' 1 1 

where we have adopted the notation (— ) x = (— )^i Xi . (As usual, the lattice spac- 
ing is taken to be unit.) An ordinary spin SU(2)s transformation corresponds to 
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right multiplication of \I/ by an SU(2)s matrix U: ^ — > \I/£/, while the pseudospin 
SU(2)c- transformation is left multiplication: \I> — > Obviously, the latter gener- 
ates the electron-hole transformation. At half-filling, defined by I] xo -n cr (x) = E x l, 
the Lagrangian can be rewritten in terms of ^: 

L = ^hr^fUdt^+t ]T tr[^t( x ,t)^(x',t)] 

x 1 (xx'J 

-^EH-r^*]) 2 ' ( 3 ) 

where a a (a = 1, 2, 3) are the Pauli matrices. The theory now is manifestly SU(2)fX 
SU(2)s invariant. The interaction term in Eq. (3) is devised in such to our later 
purpose of studying superconductivity and CDW order (although obviously it can be 
written in several equivalent forms). By Noether's method we may derive from this 
Lagrangian the conserved spin and pseudospin charges, 

S = E ^tr[*%«7] = Y \i c \ c l + h - c - , ic \ c i + h - c - , C T C T - c l c i}' ( 4 ) 
J = E ^tr[*t^] = 2 i{(-) x ( Ci c T + h.c.) , (-) x (^c iCT + h.c.) , 4 c CT - 1}, (5) 

which are the generators respectively of spin and pseudospin symmetries. (J3 may 
be identified as the ordinary V{l)c charge generator apart from a constant.) The 
quantum operators S and J satisfy the commutation relations: [S a , Sb] = ie a bcS c , 
[J a , Jb\ = itabcJc, and [S a , J b \ = 0. 

The SO (4) symmetry can be spontaneously broken in various patterns (see Ta- 
ble 1). Several authors [10,11,12] showed that spin textures (skyrmions) can exist in 
an antiferromagnetic (Neel) background on a square lattice. In fact, such skyrmions 
are associated with the spontaneously breakdown of the spin rotational symmetry 
SU(2)s to U(l)5, as exhibited by the 0(3) cr-model [13]. Yet, we shall explore an 
alternative case of SU(2) C broken to a 'hidden' U(l) subgroup (not U(l)c in general), 
for which the superconducting and CDW long-range orders coexist [5]. 

We now proceed to derive an effective field theory of low energy excitations. In 
the path-integral formalism, the action functional is given by the same form as Eq. (3) 
except replacing the electron annihilation and creation operators by the Grassmann 
fields (acting like anticommuting c- numbers): c a — > ip a and c\ — > By Hubbard- 
Stratonovich transformation we can cancel the four-fermion interaction term in Eq. (3) 
by adding to the Lagrangian a term AL = ^f^Exa(0a - |tr[(-) x ^ t a a ^]) 2 (for 
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negative U), where we have introduced three auxiliary scalar fields <p a . Then, 

+t y>i(x,t)M*,t)+h. c ] . (6) 

(XX') 

Integrating out the <p a fields shall return L to the original form. The Lagrangian (6) 
is stationary in 'field' space at the point 

0o(x,t) = \tv[{-T^^] 

= \ {Ml + h.c, itpi^ + h.c, (-) x ^Vv} . (7) 

The expression of 0o infers that is a three- vector in pseudospin space but transforms 
as a scalar under the spin group SU(2)s [14]. 

Now that the Lagrangian (6) is quadratic in fermion fields ip a and ipl, one can 
formally integrate out all of them and obtain the effective action for the fields of <j) a , 
which represent the collective modes associated with 'pseudospin' fluctuations. Along 
this direction, one can follow a standard procedure (see, e.g., Sec. 3.4 of Ref. [13]). 
However, we shall derive the low energy effective field theory with an alternative 
approach, in which the role of symmetry breaking will be evident, based on such an 
argument: the Lagrangian (6) is SU(2)c invariant, and hence any effective action 
derived from it must be also SU(2)c invariant (no matter whether realized linearly 
or nonlinearly) . Thus, together with time-reversal symmetry and parity, the SU(2)c 
invariance requires that the (continuum-limit) effective action for <p be given in the 
following general form 

1 „ ? „ ? 1 



/I I 
d d xdt - Cl $.idt$+ -c 2 d t $- dtf- -c' 2 dS- d l( f+ Y{4>) 



(8) 



r(0) = rf 2 0-0 + rf 4 (0-0) i + ---, (9) 

where space index % runs from 1 to d. The terms indicated by '• • •' will contain higher 
powers of the <fi fields and/or their derivatives. All coefficients (ci, • • •, and d 2 , ■ ■ •) 
are assumed real and finite after renormalization [15]. 

With the ground state pointing in some specific direction in pseudospin space 
(0) = A 7^ 0, the pseudospin SU(2)c symmetry is spontaneously broken down, 
whereas the ordinary spin SU(2)s rotational invariance is preserved but rather realized 
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trivially. (Recall that and A transform as scalars under SU(2)g.) By Eq. (7), 
Ai = ^^ + ^\), A 2 = l -(^-^\) and A 3 = \{{~T^), (10) 

and they are naturally interpreted as the pairing fields and CDW order parameter. 
This is a manifestation of the fact that CDW order and superconductivity can coexist 
in the negative-?/ Hubbard model at half-filling for d > 2 [5]. Nevertheless, the ground 
state of (0) = A is still invariant under SO (2) ~ U(l) rotations about the direction 

— * 

of A. This unbroken 'hidden' U(l) subgroup is however generally not the same as 
the ordinary U(l)c charge symmetry. Clearly, a transformation of SU(2) C can rotate 
the SC pairing fields A 12 and the CDW order parameter A 3 into one another [16]. 
Thus, the coexistence of SC and CDW can be perfectly understood as a result of 
the breakdown of SU(2) C — > U(l)hidden with the order-parameter space topologically 
regarded as a 2-sphere manifold S 2 . 

Since the action (8) is SU(2)c- invariant, it is always possible to rotate the pseu- 
dospin 'field' space in such a way that A aligns into the three-direction (z-direction) in 

— * 

the new 'field' space. We then express the field as a pseudospin rotation R acting on 
a three-vector (0,0, a) whose first two components vanish: a (x, t) = i? a3 (x, t)cr(x, t) 
with R ab (a, 6=1,2, 3) an orthogonal matrix -R T (x, t)i?(x, t) = 1. In place of the field 
variables <f> a , our variables now are <r(x, t) and whatever other variables are needed to 
parameterize the rotation R. Let us simply choose those parameters as the _R a3 (x, t) 
themselves but rather call them O a (x, t) hereafter. By definition, Q a represent the 
Goldstone mode degrees of freedom for the pseudospin fluctuations. Furthermore, as 
far as the lowest energy excitations are concerned, we may simply replace <r(x, t) with 
its mean-field value |A| = \/ A a A a . The effective action (8) then becomes 



I cS = J d d xdt 



C2 



with a constraint Q ■ Q = 1. This is the familiar nonlinear 0(3) cr-model. It ex- 
hibits that two gapless Goldstone modes appear around the CDW and supercon- 
ductivity background along with the spontaneous breakdown of SU(2)<7xSU(2)s — + 
U(l) xSU(2),5, as required by the Goldstone theorem [17]. We remark that the result 
is true for d > 2 space dimensions only [18]. 

The effective field theory described by Eq. (11) can predict the existence of a 
topological object, known as skyrmion. To see this, consider the energy of a static 
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solution 

E= d.\Z\ 



J d d xd i n a -d i n a + ---. (12) 
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Field configurations of finite E must have 9jQ a (x) vanishing at spatial infinity faster 
than |x|~ d / 2 (where |x| = y/XiXi), so that fi a ( x ) must approach a constant f2 aoo as 
x — > oo with a remainder vanishing faster than |x| ( 2— *0/ 2 . The fields f2 a at any 
point form a homogeneous space, the (order-parameter) coset space of SU(2)^ x 
SU(2) S /U(1) x SU(2) 5 ~ SO(3)/SO(2) = S 2 (2-sphere), for which it is possible to 
transform any one field value to any other by a transformation of SU(2)c- The field 
f2 a (x) thus represents a mapping of the whole <i-dimensional space, with the sphere 
r = oo taken as a single point, into the manifold S 2 of all field values. Therefore, finite- 
energy static configurations fi a ( x ) in d space dimensions may be classified according 
to the dth homotopy group of S 2 , ^(S^) (for a lucid discussion, see, e.g., Chap. 23 
of Ref. [17]). 

First, let us consider d — 2 case. The second homotopy group ^(S^) = Z where Z 
is the group of integers. Therefore, there exist skyrmions characterized by an integral 
winding number Q (equally called topological charge as well), which can be expressed 
as 

Q = — J d 2 xeij(l ■ (did x djQ). 

A skyrmion with winding number Q has energy E = 47ic' 2 \A\ 2 \Q\, and its analytical 
solution is already known [19]. For example, the solution for a skyrmion of winding 
number Q — 1 and scale A reads 

"i,2 (x) = , and S2 3 (x) 



x 2 + 4A 2 ' v ' x 2 + 4A 2 ' 

As of d — 3, 7r 3 (5'2) = Z. Hence, the theory also predicts the existence of pseudospin 
textures (skyrmions) classified by an integral winding number. Strictly speaking, the 
functional (12) however does not have skyrmion stationary points for d — 3, unless 
we add terms involving higher powers of dfl a to the integrand, as indicated by 
For instance, a four-derivative term will be sufficient [20]. 

One can gain an insight into the pseudospin skyrmions by comparing the pseu- 
dospin J with the ordinary spin S as defined in Eqs. (5) and (4), respectively. The 
Hubbard Lagrangian (1) is invariant under the combined action of the Lieb-Mattis 
canonical transformation defined by c-f(x) — > C|(x), Cj(x) — > (— ) x c|(x), and the re- 
versal of interaction sign U — > —U [21]. Meanwhile, the Lieb-Mattis transformation 
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— * — * — * 

maps J onto S. For U > 0, the ground state is antiferromagnetic with (S) ^ 0, 
which spontaneously breaks the ordinary spin symmetries, and the (spin) skyrmions 
are known to exist. By contrast, when U < 0, the ground state is in the phase of SC 

— * 

and CDW with (J) ^ 0, which spontaneously breaks the pseudospin symmetries, and 
consequently there appear pseudospin skyrmions. (Table 1 summaries the analogy 
between the two cases.) 

The skyrmion configuration shall become unstable if the SO (4) symmetry is ap- 
proximate. This happens when the Hubbard model is doped slightly away from 
half-filling. Then, a chemical potential term needs to be included in the Lagrangian 
(Eqs. 1 and 3): 

^brcak = ~H C l C ° = ~ 2 I J ' J 3 + Const. , (13) 
x 

which when \i ^ explicitly breaks the SO (4) symmetry down to the obvious 
U(l)cxSU(2)s subgroup. Since the Lieb-Mattis transformation maps J 3 into S3, 
this explicit breakdown of symmetry can be better viewed as caused by a 'fictitious 
magnetic field', namely, /i. The presence of L break will somehow modify the effective 
action (11) and a detailed discussion shall be given elsewhere. Nevertheless, accord- 
ing to our understanding of spontaneously broken approximate symmetries, there are 
two immediate consequences of it. One is that the Goldstone modes become massive. 
(The usual phase mode associated with U(l)c however remains massless.) Another 
is that skyrmion configurations are allowed to unwrap. Imagine to do the following 
experiment. One first creates skyrmions with some winding number in the ground 
state and then starts doping the system. We conjecture that the 'fictitious magnetic 
field' jj, should cause the skyrmions to unwrap but there will be some energy barrier 
to make it slow [22]. By contrast, the spin skyrmion texture in an antiferromag- 
netic background differs from the pseudospin skyrmion here in that the former is 
topologically stable against doping. 
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patiently teaching me homotopy theory and T.-L. Ho for stimulating conversations. 
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Table 1: Comparison of two possible patterns of symmetry breaking. Here, the Neel 
vector M = ((-) x ,S(x)) and A is defined in Eq. (10). 



Properties 


U > 


U < 


Unbroken subgroup 


SU(2) c xU(l) 5 


U(l)xSU(2) s 


Broken symmetries 


TWO of SaS 


Two of J a 's 


Phase 


Neel 


SC+CDW 


Order parameter 


\M\ 


|A| 


Skyrmion (d = 2) 


Spin texture 


Pseudospin texture 


Stability vs doping 


Stable 


Unwrap 
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